Why edge effects are important on the intrinsic loss mechanisms of graphene 

nanoresonators? 
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Molecular dynamics simulations are performed to investigate edge effects on the quality factor 
of graphene nanoresonators with different edge configurations and of various sizes. If the periodic 
boundary condition is applied, very high quality factors (3 x 10"") are obtained for all kinds of 
graphene nanoresonators. However, if the free boundary condition is applied, quality factors will 
be greatly reduced by two effects resulting from free edges: the imaginary edge vibration effect 
and the artificial effect. Imaginary edge vibrations will flip between a pair of doubly degenerate 
warping states during the mechanical oscillation of nanoresonators. The flipping process breaks 
the coherence of the mechanical oscillation of the nanoresonator, which is the dominant mechanism 
for extremely low quality factors. There is an artiflcial effect if the mechanical oscillation of the 
graphene nanoresonator is actuated according to an artificial vibration (non-natural vibration of the 
system), which slightly reduce the quality factor. The artificial effect can be eliminated by actuating 
the mechanical oscillation according to a natural vibration of the nanoresonator. Our simulations 
provide an explanation for the recent experiment, where the measured quality factor is low and 
varies between identical samples with free edges. 

PACS numbers: 62.40.+i, 63.22.Rc, 68.65.-k 
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I. INTRODUCTION 

Several recent works have investigated possible appli- 
cations of the one- atom-thick graphene in the resonator 
field. The performance of resonators is described by the 
quality factor, which reflects the energy damping of the 
mechanical oscillation for the resonator. Both optical and 
electrical methods are applied to actuate the mechanical 
oscillation of graphene nanoresonators. Applying these 
two techniques, Bunch et al. detected quite low quality 
factorSji while Zande et al. found in a recent experi- 
ment that the quality factor increases dramatically with 
cooling and can reach up to 9000 at 10 K with graphene 
resonators produced from the chemical vapor deposition 
growth method.^ A powerful technique in the nanores- 
onator field is to identify the shape of the mechanical 
oscillation of the nanoresonator by using a novel form of 
scanning probe microscopy.— This technique is useful for 
the detection of which natural vibration mode has been 
excited during the actuation of the mechanical oscillation 
of the nanoresonator. Particularly, it is useful to examine 
whether the mechanical oscillation of the nanoresonator 
is following a natural vibration of the system or not. 

The quality factor is limited by both external and in- 
trinsic damping mechanisms. Many efforts have been 
devoted to the investigation of various external damping 
mechanisms, such as charges in the substrate, attach- 
ment losses, and etc.* The intrinsic nonlinear effect was 
studied by a continuum elastic model derived from atom- 
istic interaction for the graphene nanoresonator,^ or by 
molecular dynamics simulations.^ As another important 
intrinsic energy loss mechanism, the edge effect from free 
edges in the graphene was investigated by Kim and Park 



using molecular dynamics simulations in 2009^. They 
show that edge vibrations on a pair of free edges in the 
graphene accelerate the loss of the coherence of the me- 
chanical oscillation of the nanoresonator, which results 
in very low quality factors. The recent experiment by 
Zande et al. also found that the edge effect is impor- 
tant for the variations in the resonance frequency and 
the quality factor of the graphene nanoresonator)^ which 
has not received adequate interpretation. This experi- 
ment has also found that free edges in the graphene can 
greatly reduce quality factors. 

In recent years, much attention has been paid to the 
configuration of free edges in graphene nanoribbon. Gass 
et al. apply the high-angle annular dark-field tech- 
nique to observe the edge reconstruction in graphene 
nanoribbon^ An unusual scrolling and staggering process 
was observed during the reconstruction. This unusual 
process is eventually attributed to the warping states, 
which are induced by intrinsic edge stress in graphene.— 
Warping states are localized edge states with out-of-plane 
movement, and they are also found to be responsible for 
the edge reconstruction of the graphene nanoribbon un- 
der Joule heating treatment J^iii Now, a natural question 
arises: is there any relation between these warping states 
and the reduction of quality factors in graphene nanores- 
onators by free edges? 

In this paper, we report molecular dynamics simula- 
tions to study edge effects on graphene nanoresonators 
with both zigzag and armchair edges, and of various 
sizes. Very high quality factors are obtained for graphene 
nanoresonators with periodic boundary condition (PBC), 
and the quality factor is inversely proportional to temper- 
ature as 1/T, which is due to phonon-phonon scattering. 
However, quality factors have quite different value for 
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FIG. 1: (Color online) Energy time history of Sample- A at 
10 K with PBC for two short edges, (a) Kinetic energy, (b) 
Potential energy, (c) Close-up of kinetic energy, (d) Close-up 
of potential energy. 

zigzag and armchair graphene nanoresonators with free 
boundary condition (FBC). For zigzag edges, there are 
imaginary edge vibrations localized at free edges, which 
will flip between two doubly degenerate warping states 
during the mechanical oscillation of graphene nanores- 
onators. This flipping process breaks the coherence of 
the mechanical oscillation of the graphene nanoresonator 
and leads to extremely low quality factors with tempera- 
ture dependence 1/T°-^®. The quality factor of zigzag 
graphene nanoresonator will be further reduced with 
temperature dependence if the oscillation of the 

resonator is actuated according to an artificial vibra- 
tion (non-natural vibration of the system) . For armchair 
edges, there is no imaginary edge vibration, so quality 
factors are mainly reduced by the artificial effect. An ex- 
tremely high quality factor can be obtained for armchair 
grphene nanoresonators if the mechanical oscillation is 
actuated naturally (according to its natural vibration), 
although FBC is applied. 

The rest of the paper is organized as follows. Sec. 11(A) 
is devoted to some techniques for the simulation and cal- 
culation approach. In Sec. 11(B), the imaginary edge vi- 
bration effect on quality factors is discussed. The arti- 
ficial effect on quality factors is illustrated in Sec.II(C). 
The paper is concluded in Sec. III. 

II. RESULTS AND DISCUSSION 

A. simulation setup 

The first simulated graphene sample is of size 
{L^,Ly) = (19.7,127.8) A, which will be referred to 
as Sample- A throughout the text. Short edges are two 
zigzag edges, so this type of graphene is usually called 
zigzag graphene nanoribbon. The x-axis is along the 
short edge, and the y-axis is along the long edge. The 
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FIG. 2: (Color online) Kinetic energy time history of Sample- 
A at 10 K with FBC. The mechanical oscillation is actuated 
artificially in (a) and naturally in (b). 



z-axis is perpendicular to the graphene plane. Two long 
edges are always fixed in all simulations. We apply either 
PBC or FBC for the two short edges. The interatomic 
interaction is described by the Brenner potential.''^ The 
Newton equations of motion are integrated using the ve- 
locity Verlet algorithm with a time step of 1 fs. 

Fig. [1] shows the time history of kinetic energy in panel 
(a) and potential energy in panel (b) for Sample-A at 10 
K. PBC is applied at the two short edges. For i < 0, 
the Nose-Hoover-i^iii heat bath is applied to thermalize 
the system to a constant temperature within the NVT en- 
semble. In panel (a), there is a 5% variation in the kinetic 
energy due to the statistical natural of the Nose-Hoover 
algorithm. At t = 0, the nanoresonator is actuated by 
adding a sinuous velocity Avj oc Cz sin^nxj / L^) to each 
atom j, where Xj is the a;-coordinate of the atom and Cz is 
the unit vector in z direction. The energy increase by the 
added velocity Av is 50% of the total kinetic energy. It is 
large enough to ensure that the mechanical oscillation of 
the nanoresonator is obviously stronger than the thermal 
vibration in the nanoresonator. Yet, it is less than 0.003% 
of the total potential energy. For t > 0, the nanores- 
onator oscillates within the NVE ensemble. Panel (c) / (d) 
is the close-up of the kinetic/potential energy. Panels (c) 
and (d) demonstrate the energy exchange between kinetic 
and potential energy after t > 0, while the total energy is 
kept as a constant. The energy exchange between kinetic 
and potential energy reflects the mechanical oscillation of 
the nanoresonator. Some roughness can be observed on 
top of the curve in panels (c) and (d), which display the 
thermal vibration at 10 K. The decay of the oscillation 
amplitude is used to analysis the quality factor of the 
nanoresonator by fitting the kinetic energy after i > to 
a function (t) = a-|-5(l — 27r/(3)*cos(a;t). w is the nat- 
ural frequency of the nanoresonator. Parameter a gives 
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FIG. 3: (Color online) Two edge vibrations localized nearby 
the two free short edges in Sample-A. The arrow attached 
on each atom respresents the vibrational amplitude of this 
atom in the vibration. Vibrational amplitudes decay expo- 
nentially from edge into center. The frequencies of these two 
edge vibrations in (a) and (b) are degenerate and imaginary, 
= -625 cm~^ 



final kinetic energy after the mechanical oscillation dies 
away. The decay of the mechanical oscillation will lead 
to the temperature increase in the nanoresonator. Q is 
the resulted quality factor. The symbol Q will be used 
to denote the normal mode coordinate throughout the 
text. This is the only place where Q denotes the quality 
factor. The obtained quality factor is as high as 3 x 10^ 
at 10 K. Similar fitting procedure can be conducted for 
the potential energy, which results in a same value for the 
quality factor. In the following, the quality factor will be 
calculated from the fitting of the kinetic energy. 



B. effects from imaginary edge vibrations 

It was shown that the quality factor of a graphene 
nanoresonator will be dramatically reduced by free 
edges. ^ We repeat this result in Fig. [5] (a) for Sample-A, 
which is also the sample studied by Kim and Park. FBC 
has been applied for the simulation, while all other pa- 
rameters are the same as Fig. [1] In particular, we would 
like to stress that the mechanical oscillation is also actu- 
ated by adding the sinuous velocity Aw, which is the same 
as Ref.0. Our simulation demonstrates that the mechan- 
ical oscillation decays very fast in presence of free edges, 
resulting in an extremely low quality factor. 

The origin for this edge effect is still unclear, al- 
though the phenomenon of low quality factor due to 
free edges has been shown in previous work. To reveal 
the underlying mechanism, we study all natural vibra- 
tions in Sample-A with FBC by solving the eigenvalue 
problem of dynamical matrix derived from the Bren- 
ner potential. The dynamical matrix is obtained from 
Kij = d^V/dxidxj, where V is the Brenner potential 
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FIG. 4: Brenner potential energy for graphene nanoribbon 
with position fj ~ + ^Xcdgcvibration ^'^^ atom j, where is 
the optimized structure. Xcdgcvibration is the eigen vector of 
the imaginary edge vibration localized at one short edge as 
shown in Fig. [3] (a). The right top inset displays a configu- 
ration for > 0, where edge atoms are above the graphene 
plane. The left bottom inset demonstrates a configuration for 
77 < 0, where edge atoms are below the graphene plane. 



and Xi is the position of the i-th degree of freedom. This 
formula is realized numerically by calculating the energy 
change after displace the i-th and j-th degrees of freedom 
for a small value. We find several edge vibrations local- 
izing at the two free edges. Fig. [3] displays two lowest- 
energy edge vibrations. The arrow attached on each atom 
represents the vibrational amplitude of this atom in the 
vibration. Vibrational amplitudes decay exponentially 
from edge into center, exhibiting the localization prop- 
erty. These two edge vibrations are degenerate. They 
localize at each end of the two free edges. An interest- 
ing property of these two edge vibrations is their imag- 
inary frequency, = —625 cm~^. We have optimized 
the structure before the derivation of the dynamical ma- 
trix from the Brenner potential. During the structure 
optimization, the residue force on each atom has been 
required to be as small as 10"'* eV/A, so these two imag- 
inary vibrations do not result from numerical errors; in- 
stead, they refiect the intrinsic properties of free zigzag 
edges. 

Imaginary edge vibrations imply the instability of free 
zigzag edges when the graphene nanoribbon vibrates ac- 
cording to these two imaginary edge vibrations. Simi- 
lar instability was also found by Lan et al. using tight- 
binding potentialJ^ Starting with these imaginary vibra- 
tions, we can search for the ground state correspond- 
ing to a stable structure. Fig. |3] shows the potential 
energy of the graphene nanoribbon with position r^- = 

^ + ^Xodgovibration ^r atom j, where r° is the optimized 
structure. Xodgovibration is the eigen vector of the imagi- 
nary edge vibration localized at one short edge as shown 
in Fig. [3] (a). Similar potential curve is also found for 
the other imaginary edge vibration shown in Fig. |3] (b). 
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FIG. 5: (Color online) Normal mode coordinate time history 
of the two edge vibrations in Fig. [3] showing flipping phenom- 
ena. The simulated system is Sample-A at 10 K with FBC. 
The mechanical oscillation is actuated artificially, (a) is for 
the edge vibration shown in Fig. |3] (a), (b) is for the edge 
vibration shown in Fig. [31(b). (c) is the close-up for (a), (d) 
is the close-up for (b). 



The parameter 77 is in [-1.5, 1.5] in the figure. Two insets 
display the edge configuration for the graphene nanorib- 
bon with 77 > (right top inset) and 77 < (left bot- 
tom inset). The potential curve demonstrates that the 
optimized structure (at 77 = 0) actually corresponds to 
a local maximum value of potential energy instead of a 
minimum value. Meanwhile, there are two ground states 
at 77 = ±0.9, which are doubly degenerate to each other. 
We refer to these two ground states as 771 and 774, states, 
respectively. They are the warping states, where the 
graphene nanoribbon is scrolled and warped at the edge.^ 
The graphene nanoribbon has a mirror reflection symme- 
try, Gxy, with respect to the xy plane, resulting from its 
two-dimensional nature. Each of the warping state alone 
breaks this mirror reflection symmetry as can be seen 
from both insets. This is a spontaneous symmetry break- 
ing phenomenon resulting from the Mexican-top-hat-like 
potential shown in Fig. 21— i.e one symmetrical system 
ends up with an apparently asymmetric ground state, 
which is ubiquitous in condense matter physics such as 
superconductorsii and superfluidsr^, and is related to 
the Nambu-Goldstone boson in the scalar field theory^^. 
The mirror reflection symmetry can be restored if these 
two warping states are considered together, as the mir- 
ror operation transforms one warping state to the other. 
The optimization at K preserves the mirror reflection 
symmetry in the graphene nanoribbon; thus it can not 
achieve either of the two ground (warping) states. The 
frequency of the imaginary vibration (25i cm"-'^) gives us 
some information about the energy barrier between these 
two warping states: fioj — 25/8.06554 « 3.1 meV. This 
value is comparable with the energy barrier, 5.0 meV, as 
reading directly from the potential curve. 

As the edge vibration manifests the characteristic dif- 
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FIG. 6: The edge configuration of the graphene nanoresonator 
undergoes a sudden change from r]i to 77^ state during the first 
flipping process in Fig. [5](c). 

ference between PBC and FBC, it is natural to attribute 
the extremely low quality factor to edge vibrations. To 
verify this thinking, we monitor the evolution of these 
two imaginary edge vibrations during the molecular dy- 
namics simulation. The normal mode coordinate for a 
natural vibration k can be obtained from its standard 
definition 

Qk{t)^Y.XjirAt)~r'^), (1) 

where 3A'' is the total degrees of freedom, x*^ is the eigen 
vector for this natural vibration. rj{t) is the position 
of the j-th degree of freedom and is the equilibrium 
position. Time histories of Q{t) for these two imaginary 
edge vibrations are shown in Fig. [5] For t < 0, the ther- 
mal process is controlled by heat bath and Q{t) oscillates 
around a nonzero value, which is a more stable state. Af- 
ter the removal of heat bath for t > 0, there is a flipping 
phenomenon; i.e the normal mode suddenly flips between 
two stable states. The magnitudes of Q{t) are the same 
in these two stable states, while the signs are opposite. It 
indicates that energy does not change after the flipping 
process, since energy is proportional to Q^{t). We now 
justify that these two stable states are exactly the two 
warping states discussed above. The two insets in Fig. [S] 
show the edge configurations of the graphene nanores- 
onator before and after the flipping process. Indeed, the 
system jumps from the warping state 7/4^ to the other 
warping state rj^ during this flipping process. These two 
warping states are degenerate, so the energy does not 
change during the flipping process. This corresponds to 
the above fact that the Q{t) changes sign with unchanged 
amplitude during the flipping process. The flipping pro- 
cess only exists in the early time history, when the system 
is out of thermal equilibrium. There is no flipping process 
after the mechanical oscillation dies away; i.e the whole 
system reaches thermal equilibrium again. 
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FIG. 7: (Color online) Kinetic energy time history at 300 K 
for Sample-A with (a) PBC, (b) FBC and artificial actuation, 
and (c) FBC and natural actuation. 




FIG. 8: (Color online) Quality factor versus temperature for 
Sample-A with PBC (black filled circle), FBC and artificial 
actuation (red filled square), and FBC and natural actuation 
(blue filled trigon). 



The flipping phenomenon is a disaster for the quality 
factor of the mechanical resonator, because the positions 
of atoms in the edge region are suddenly relocated during 
each flipping. This is extremely harmful for the coherence 
of the mechanical oscillation, leading to very low quality 
factors. Panels (a) and (b) in Fig. [7] are simulation re- 
sults for Sample-A at 300 K, which are similar to that 
at 10 K. The energy damping for the nanoresonator with 
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FIG. 9: (Color online) Kinetic energy time history at 10 
K for zigzag graphene nanoresonator of size {Lx,Ly) — 
(19.7, 29.82) A with (a) PBC, (b) FBC and artificial actu- 
ation, and (c) FBC and natural actuation. 



FBC is obviously much stronger than that with PBC. 
Fig. [5] shows the temperature dependence for the quality 
factor. It shows that the quality factor for the nanores- 
onator with FBC (red flUed squares) is about three orders 
smaller than the nanoresonator with PBC (black filled 
circles) in whole temperature range. The quality factor 
for the nanoresonator with PBC is inversely proportional 
to the temperature as 1/T. This temperature behav- 
ior results from the three-phonon scattering, which is 
the only intrinsic decaying mechanism in the graphene 
nanoresonator with PBC. For nanoresonators with FBC, 
the power factor in the temperature dependence is con- 
siderably reduced to be only 0.16. It is because the imag- 
inary edge vibration effect (flipping process) is not sensi- 
tive to the temperature. The combination of the imagi- 
nary edge vibration effect and phonon-phonon scattering 
effect results in a power factor smaller than 1. A small 
power factor of 0.16 indicates that the imaginary edge 
vibration effect overcomes the phonon-phonon scattering 
effect. 

Fig. ini shows time history of kinetic energy for zigzag 
graphene nanoresonator of another size {Lx,Ly) = 
(19.7, 29.82) A at 10 K. There are also two imaginary 
edge vibrations at the free zigzag edges of this graphene 
nanoribbon. The frequencies of these imaginary edge vi- 
brations are almost the same as that of Sample-A, be- 
cause edge vibrations are only sensitive to the configura- 
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FIG. 10: (Color online) Two edge vibrations localized nearby 
the two free short edges in armchair graphene nanoresonator 
of size {Lx,Ly) = (25.56,49.2) A. Vibrational amplitudes de- 
cay exponentially from edge into center. The frequencies of 
these two edge vibrations in (a) and (b) are degenerate at 
3.02 cm"\ 
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FIG. 11: (Color online) Normal mode coordinate time history 
of the two edge vibrations in Fig. [10] for armchair graphene 
nanoresonator of size (L^, Ly) = (25.56, 49.2) A with FBC at 
10 K. (a) and (b) are for these two edge vibrations in artificial 
resonator, (c) and (d) are for natural resonator. 



tion of edge regions which are the same in both systems. 
We observe similar phenomenon in this nanoresonator. 
Panel (a) shows that there is only weak energy damping 
in the nanoresonator with PBC, resulting in a very high 
quality factor. Panel (b) shows much stronger energy 
damping for the nanoresonator with FBC, which leads 
to an extremely low quality factor. 

In the above, we find imaginary edge vibrations in 
zigzag graphene nanoresonators with FBC. Our simula- 
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FIG. 12: (Color online) Kinetic energy time history at 10 K 
for armchair graphene resonator of {Lx,Ly) = (25.56,49.2) A 
with (a) PBC, (b) FBC and artificial actuation, and (c) FBC 
and natural actuation. 



tions have disclosed a fact that the flipping processes of 
these imaginary edge vibrations break the coherence of 
the resonator's mechanical oscillation, which leads to ex- 
tremely low quality factors in zigzag graphene nanores- 
onators with FBC. It is natural to ask what will be the 
quality factor for graphene nanoresonators without imag- 
inary edge vibration. If our above analysis based on the 
flipping processes of imaginary edge vibrations are true, 
we should obtain weak energy damping and high qual- 
ity factors for resonators without imaginary edge vibra- 
tion. To confirm this assumption, we carry out simula- 
tion for graphene nanoresonators with armchair edges, 
where there is no imaginary edge vibration. The first 
armchair graphene nanoresonator we studied is of size 
{L^,Ly) = (25.56, 49.2) A. Fig.HUldisplay the two lowest- 
energy edge vibrations for FBC. The vibrational ampli- 
tude also shows an exponential decaying from edge into 
center region, which is a characteristic of localized vi- 
brations. These two edge vibrations are degenerate at 
3.02 cm~^ because of two free edges, and their eigen 
vectors happen to mix together. Normal mode coordi- 
nates of these two edge vibrations at 10 K are demon- 
strated in Fig. [TT] (a) and (b). Obviously, there is no 
flipping phenomenon anymore. Instead, these two edge 
vibrations involve in continuous vibrations, which helps 
to preserve the coherence of the resonator's mechanical 
oscillation. As a result, we should observe very weak en- 
ergy damping. Indeed, Fig. [T^] shows that the energy 
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FIG. 13: (Color online) Kinetic energy time history at 
10 K for armchair graphene nanoresonator of {Lx,Ly) — 
(25.56, 24.6) A with (a) PBC, (b) FBC and artificial actu- 
ation, and (c) FBC and natural actuation. 



damping for the mechanical oscillation is very weak in 
resonators with both PBC (in panel (a)) and FBC (in 
panel (b)). It means that the quality factor is not nec- 
essarily low for graphene nanoresonators even if FBC is 
applied. These results provide an indirect evidence for 
the correctness of our analysis of attributing low quality 
factor to the flipping process of imaginary edge vibra- 
tions. Fig. [13] shows similar kinetic energy time history 
for nanoresonator with armchair edges of another size 
{{L^,Ly) = (25.56,24.6) A) at 10 K. 



C. effects from artificial or natural actuation 

Up to now, we have revealed one underlying mecha- 
nism for extremely low quality factors in zigzag graphene 
nanoresonators with free edges, where the flipping pro- 
cesses of imaginary edge vibrations play an important 
role. We have also shown that high quality factors can 
be obtained for armchair graphene nanoresonators with 
free edges, because of the absence of imaginary edge vi- 
brations. Now we switch to the discussion of another ef- 
fect from free edges in graphene nanoresonators. From a 
basic point of view, the mechanical oscillation of a natural 
resonator is actuated by mechanically exciting a partic- 
ular natural vibration of the system. After the actuation 
of the resonator, the frequency of the resonator equals 
to the frequency of the natural vibration, and the me- 
chanical oscillation of the resonator corresponds to the 




FIG. 14: (Color online) (a) and (b) are two different side 
views for the first natural vibration (non-edge vibration) of 
Sample-A with PBC. It is only sinuous function in x direction, 
(c) and (d) are two different side views for the first natural 
vibration (non-edge vibration) of Sample- A with FBC. It is 
sinuous function in both x and y directions. 
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FIG. 15: The fitting of the first natural vibrations to sinu- 
soidal function y — a * sm(x — h)/c for Sample-A. (a), is 
for the distribution of the vibration in x direction with pe- 
riodic boundary condition applied in y direction, resulting 
in a = 0.053, h = 2.64 A, and c = 4.58 A. (b). is for the 
distribution of the vibration in x direction with free bound- 
ary condition applied in y direction, resulting in a = 0.082, 
h = 2.64 A, and c = 4.58 A. (c). is for the distribution of 
the vibration in y direction with free boundary condition ap- 
plied in y direction, resulting in a = 0.078, b — 10.35 A, and 
c = 38.91 A. 



vibrational morphology of the natural vibration. Let's 
review how we actuated nanoresonators in the above. In 
all previous simulations, nanoresonators are actuated by 
adding a velocity Avj oc ezSin{xj/Lx) to each atom j. 
Fig. [M] (a) and (b) show that this sinuous function is a 
natural vibration in Sample-A with PBC. We use (with 
j — 1, 2, 3, 3iV) to denote eigen vectors of all natural 
vibrations in graphene nanoresonators with PBC; partic- 
ularly, we use — sin(x/La;). We use a sinusoidal func- 
tion y = a shi{x — b) j c to fit mode in Fig. [TS] (a) , where 
the relative fitting errors for parameters (a, h and c) are 
less than 0.1%. This vibration does not depend on the y 
coordinate, so all data falls in one sinusoidal function. It 
should be noted that the sinuous function is the natural 
vibration in a rectangular plate, while the Bessel function 
is the natural vibration in the circular plate (such as the 
circular graphene in Ref. [2ll))^ The graphene nanores- 
onator with PBC is a natural resonator if its mechani- 
cal oscillation is actuated according to ^i. However, the 
sinuous function is not a natural vibration in Sample-A 
with FBC. We use Xj (with j = 1, 2, 3, 37V) to denote 
eigen vectors of all natural vibrations in Sample-A with 
FBC. The lowest-energy natural vibration (non-edge vi- 
bration) is shown in Fig. [14] (c) and (d) for Sample-A with 
FBC. We denote this particular natural vibration as xi- 
Fig. [TH] (b) and (c) are the fitting for the distribution of 
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FIG. 16: (Color online) (a) Decomposition of the first natu- 
ral vibration of Sample- A with PBC (^i) in terms of the nat- 
ural vibrations of Sample-A with FBC (xj, j=l,2,3,...,3Af): 
C'lj = (?i|Xj)- Two blue solid arrows indicate the two 
edge vibrations shown in Fig. |3] for Sample-A with FBC. 
The imaginary value of the frequency uj = ±25icm~^ has 
been intensively shown as a negative value to = — 25cm~^. 
The red dotted arrow indicates the first natural vibration of 
Sample-A with FBC (xi)- (b) Similar decomposition for the 
first natural vibration of armchair graphene nanoresonator of 
{L^,Ly) = (25.56,49.2) A with PBC. Two blue solid arrows 
indicate the two edge vibrations shown in Fig. [TOl 



(a) 













1 2 





(b) 1 








1 2 3 0.2 0.4 

time (ns) time (ns) 



FIG. 17: (Color onfine) Normal mode coordinate time history 
of the two edge vibrations in Fig. [3] showing flipping phenom- 
ena. The simulated system is Sample-A at 10 K with FBC 
and natural actuation, (a) is for the edge vibration shown in 
Fig. [3] (a), (b) is for the edge vibration shown in Fig.|3](b). (c) 
Close-up for the flipping phenomenon in the flrst edge vibra- 
tion, (d) Close-up for the flipping phenomenon in the second 
edge vibration. 



vibration of xi- The fitting curve in panel (b) is for the 
those atoms in the middle of the graphene nanoribbon 
vifith y « 70 A, while the fitting curve in panel (c) is for 
the those atoms in the middle of the graphene nanorib- 
bon with X « 10 A. It clearly shows that this natural 
vibration is a sinuous function in both x and y direc- 
tions. Hence, the graphene nanoresonator with FBC is 
not a natural resonator if its mechanical oscillation is ac- 
tuated according to ^i. We name this type of resonator 
as an artificial resonator, which is actuated according to 
an unnatural vibration of the system. For an artificial 
resonator, many natural vibrations are excited simulta- 
neously during the actuation of the mechanical oscillation 
of the nanoresonator. can be decomposed in terms of 
Xj (with J = 1,2,3,..., 3iV): 



3N 



(2) 



Those vibrations with large coefficients \cij\ are consid- 
erably excited during the actuation of the nanoresonator 
according to ^i. Fig. [12] (a) displays the coefficient cij 
for all natural vibrations in Sample-A with FBC. Those 
natural vibrations above 100 cm~^ are not excited. The 
highest peak (indicated by red doted arrow) corresponds 
to the natural vibration xi in Sample-A with FBC. We 
find that the other two strong peaks (indicated by blue 
solid arrows) are with respect to the two imaginary edge 
vibrations at the free edges as shown Fig. [31 Actually, 
from the vibrational morphology of natural vibrations , 
Xi , and the two imaginary edge vibrations 

Xcdgcvibration ; 

the relationship among them can be simply obtained: 
~ Xi + Xcdgcvibration- This dccomposition indicates 
that the two edge vibrations are also considerably excited 
when the nanoresonator with FBC is actuated according 
to ^1. These edge vibrations are not in thermal equi- 
librium state. As a result, the edge effect is stronger in 
artificial resonator than the natural resonator. 

The imaginary edge vibration effect (discussed in pre- 
vious section) and the artificial effect are both important 
mechanisms for low quality factors in graphene nanores- 
onators with FBC. We can split these two effects sep- 
arately. There is no artificial effect if the graphene 
nanoresonator with FBC is actuated naturally accord- 
ing to xi- This is a natural mechanical resonator, so it 
can be used to study the pure imaginary edge vibration 
effect. Fig. [2] (b) shows the time history of kinetic en- 
ergy in Sample-A with FBC, which is naturally actuated 
according to its xi- The energy damping is still very 
strong, which is only slightly smaller than the artificial 
resonator shown in panel (a) of the same figure. The nor- 
mal mode coordinates for the two imaginary edge vibra- 
tions in Samplc-A arc monitored in Fig. 1171 We observe 
the same flipping phenomenon as that in the artificial 
resonator shown in Fig. [5] As we have discussed above, 
the flipping process breaks the coherence of the mechani- 
cal oscillation, leading to a low quality factor. The strong 
energy damping in Fig. [2| (b) means that the imaginary 
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edge vibration effect is the key mechanism for low qual- 
ity factors; while the small difference between Fig. [5] (b) 
and (a) indicates that the artificial effect is only of minor 
importance for the reduction of quality factor. The qual- 
ity factor is mainly determined by the imaginary edge 
vibration effect in graphene nanoresonators with imagi- 
nary edge vibrations. Fig.[5]shows that the quality factor 
in the natural resonator is slightly larger than the artifi- 
cial resonator. The temperature dependence power fac- 
tor in the natural resonator is 0.28, which is larger than 
the value of 0.16 for the artificial resonator. This slight 
difference is due to the artificial effect. The kinetic en- 
ergy time history for zigzag graphene natural resonator 
with FBC of another size is shown in Fig. |9] (c), which 
confirms that the imaginary edge vibration effect is the 
major reason for low quality factors. 

We have previously shown that there is no imaginary 
edge vibration in the armchair graphene nanoresonator 
with FBC, so it can be used to study the pure artificial 
effect on the quality factor. Fig. [T^] (c) shows the kinetic 
energy time history for armchair graphene nanoresonator 
with FBC of size {L^,Ly) = (25.56, 49.2) A. The mechan- 
ical oscillation of the nanoresonator is actuated naturally 
according to xi. The result displays very high quality 
factor in this natural mechanical resonator, although the 
FBC has been applied. The minimum energy damping is 
as weak as the mechanical resonator with FBC shown in 
panel (a) . Fig. [TT] (c) and (d) monitor the time history 
of normal mode coordinates of the two edge vibrations in 
this system. First of all, there is no hipping phenomenon 
due to the absence of imaginary edge vibrations. Fur- 
thermore, Q{t) is further suppressed after the heat bath 
is removed at t — 0, which indicates that these two edge 
vibrations are quite difficult to be excited if one does not 
excite them artificially. Fig. [TT] (a) and (b) show that 
they will be considerably excited if the mechanical os- 
cillation of the graphene nanoresonator is acuated artifi- 
cially according to which can be decomposed in terms 
of Xk in Fig. [1^] (b). Fig. [T3] (c) shows same results in 
armchair graphene nanoresonator with FBC of another 
size {Lx,Ly) = (25.56,24.6) A, which is also actuated 
naturally. 

Based on the two edge effects discussed above, we 
can give some explanations for the experiment. In the 
experimental it was found that edge effects are impor- 
tant for quality factors of graphene nanoresonators. For 
samples with free edges, the measured quality factors 
are low and the values vary between identically prepared 
samples. We can learn from our simulations that the 
imaginary edge vibrations are of key importance for low 
quality factors. Another experiment has shown that it is 
easily for free zigzag edges in the graphene nanoribbon 
to undergo a reconstruction and become more stablci^iiS 
The reconstruction is attributed to the edge vibrations 



at the free edge^^ thus it will partially eliminate the 
imaginary edge vibrations. Considering the later exper- 
iment, a reasonable explanation for the experiment in 
Ref. 2 is that some samples happen to reconstruct their 
free edges, which help to remove some imaginary edge vi- 
brations. As a result, the measured quality factor will be 
higher in these samples undergoing reconstruction, while 
quality factors are still low in samples without recon- 
struction. This explains the variation of quality factors 
between identically prepared graphene nanoresonators. 



III. CONCLUSION 

Our molecular dynamics simulations have disclosed 
two effects from the free edges on the graphene nanores- 
onator, which are both responsible for low quality fac- 
tors. They are the imaginary edge vibration effect and 
the artificial effect. Imaginary edge vibrations exist at 
free zigzag edges, so the imaginary edge vibration ef- 
fect is of crucial importance for low quality factors in 
zigzag graphene nanoresonators; while the artificial effect 
only slightly further reduces quality factors. The hipping 
process for imaginary edge vibrations between two dou- 
bly degenerate warping states will seriously break the 
coherence of the mechanical oscillation of the graphene 
nanoresonator. The power factor in the temperature de- 
pendence of quality factors is 1.0 for graphene nanores- 
onator with FBC, which reffects the intrinsic phonon- 
phonon scattering mechanism for the energy damping. 
This power factor decreases to be 0.28 in natural nanores- 
onators with FBC. This reduction is due to the hipping 
phenomenon for the imaginary edge vibrations. It is fur- 
ther reduced to be 0.16 in artificial nanoresonators with 
FBC. This further reduction results from the artificial 
effect. 

There is no imaginary edge vibration in armchair 
graphene nanoresonator with FBC, so the artificial effect 
is the only mechanism for the decrease of quality factors. 
An extremely high quality factor can be obtained for arm- 
chair graphene nanoresonator with FBC, if its mechani- 
cal oscillation is naturally actuated. This result indicates 
that free edges not necessarily lead to low quality factors 
for graphene nanoresonators. Our simulations provide 
advices for experimentalists that high quality factors can 
be obtained if graphene nanoresonators with FBC are of 
armchair edges and actuated naturally. 
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